Search for: All records

Context.Mapping molecular line emission beyond the bright low-JCO transitions is still challenging in extragalactic studies, even with the latest generation of (sub-)millimetre interferometers, such as ALMA and NOEMA.

Aims.We summarise and test a spectral stacking method that has been used in the literature to recover low-intensity molecular line emission, such as HCN(1−0), HCO⁺(1−0), and even fainter lines in external galaxies. The goal is to study the capabilities and limitations of the stacking technique when applied to imaged interferometric observations.

Methods.The core idea of spectral stacking is to align spectra of the low S/N spectral lines to a known velocity field calculated from a higher S/N line expected to share the kinematics of the fainter line (e.g. CO(1−0) or 21 cm emission). Then these aligned spectra can be coherently averaged to produce potentially high S/N spectral stacks. Here we used imaged simulated interferometric and total power observations at different S/N levels, based on real CO observations.

Results.For the combined interferometric and total power data, we find that the spectral stacking technique is capable of recovering the integrated intensities even at low S/N levels across most of the region where the high S/N prior is detected. However, when stacking interferometer-only data for low S/N emission, the stacks can miss up to 50% of the emission from the fainter line.

Conclusions.A key result of this analysis is that the spectral stacking method is able to recover the true mean line intensities in low S/N cubes and to accurately measure the statistical significance of the recovered lines. To facilitate the application of this technique we provide a public Python package, called PYSTACKER.

 

Carbon monoxide (CO) emission constitutes the most widely used tracer of the bulk molecular gas in the interstellar medium (ISM) in extragalactic studies. The CO-to-H 2 conversion factor, α 12 CO(1−0) , links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of α 12 CO(1−0) across all environments within and across galaxies as a function of metallicity, molecular gas opacity, line excitation, and other factors. Using spectral line observations of CO and its isotopologues mapped across a nearby galaxy, we can constrain the molecular gas conditions and link them to a variation in α 12 CO(1−0) . Here, we present new, wide-field (10 × 10 arcmin 2 ) IRAM 30-m telescope 1 mm and 3 mm line observations of 12 CO, 13 CO, and C 18 O across the nearby, grand-design, spiral galaxy M101. From the CO isotopologue line ratio analysis alone, we find that selective nucleosynthesis and changes in the opacity are the main drivers of the variation in the line emission across the galaxy. In a further analysis step, we estimated α 12 CO(1−0) using different approaches, including (i) via the dust mass surface density derived from far-IR emission as an independent tracer of the total gas surface density and (ii) local thermal equilibrium (LTE) based measurements using the optically thin 13 CO(1–0) intensity. We find an average value of ⟨ α 12 CO(1 − 0) ⟩ = 4.4  ±  0.9  M ⊙  pc −2  (K km s −1 ) −1 across the disk of the galaxy, with a decrease by a factor of 10 toward the 2 kpc central region. In contrast, we find LTE-based α 12 CO(1−0) values are lower by a factor of 2–3 across the disk relative to the dust-based result. Accounting for α 12 CO(1−0) variations, we found significantly reduced molecular gas depletion time by a factor 10 in the galaxy’s center. In conclusion, our result suggests implications for commonly derived scaling relations, such as an underestimation of the slope of the Kennicutt Schmidt law, if α 12 CO(1−0) variations are not accounted for. 

Determining how the galactic environment, especially the high gas densities and complex dynamics in bar-fed galaxy centers, alters the star formation efficiency (SFE) of molecular gas is critical to understanding galaxy evolution. However, these same physical or dynamical effects also alter the emissivity properties of CO, leading to variations in the CO-to-H₂conversion factor (α_CO) that impact the assessment of the gas column densities and thus of the SFE. To address such issues, we investigate the dependence ofα_COon the local CO velocity dispersion at 150 pc scales using a new set of dust-basedα_COmeasurements and propose a newα_COprescription that accounts for CO emissivity variations across galaxies. Based on this prescription, we estimate the SFE in a sample of 65 galaxies from the PHANGS–Atacama Large Millimeter/submillimeter Array survey. We find increasing SFE toward high-surface-density regions like galaxy centers, while using a constant or metallicity-basedα_COresults in a more homogeneous SFE throughout the centers and disks. Our prescription further reveals a mean molecular gas depletion time of 700 Myr in the centers of barred galaxies, which is overall three to four times shorter than in nonbarred galaxy centers or the disks. Across the galaxy disks, the depletion time is consistently around 2–3 Gyr, regardless of the choice ofα_COprescription. All together, our results suggest that the high level of star formation activity in barred centers is not simply due to an increased amount of molecular gas, but also to an enhanced SFE compared to nonbarred centers or disk regions.

 

We use new HCN(1–0) data from the ACA Large-sample Mapping Of Nearby galaxies in Dense gas (ALMOND) survey to trace the kpc-scale molecular gas density structure and CO(2–1) data from the Physics at High Angular resolution in Nearby GalaxieS–Atacama Large Millimeter/submillimeter Array (PHANGS–ALMA) to trace the bulk molecular gas across 25 nearby star-forming galaxies. At 2.1 kpc scale, we measure the density-sensitive HCN/CO line ratio and the star formation rate (SFR)/HCN ratio to trace the star formation efficiency in the denser molecular medium. At 150 pc scale, we measure structural and dynamical properties of the molecular gas via CO(2–1) line emission, which is linked to the lower resolution data using an intensity-weighted averaging method. We find positive correlations (negative) of HCN/CO (SFR/HCN) with the surface density, the velocity dispersion, and the internal turbulent pressure of the molecular gas. These observed correlations agree with expected trends from turbulent models of star formation, which consider a single free-fall time gravitational collapse. Our results show that the kpc-scale HCN/CO line ratio is a powerful tool to trace the 150 pc scale average density distribution of the molecular clouds. Lastly, we find systematic variations of the SFR/HCN ratio with cloud-scale molecular gas properties, which are incompatible with a universal star formation efficiency. Overall, these findings show that mean molecular gas density, molecular cloud properties, and star formation are closely linked in a coherent way, and observations of density-sensitive molecular gas tracers are a useful tool to analyse these variations, linking molecular gas physics to stellar output across galaxy discs.

 

Abstract We compare mid-infrared (mid-IR), extinction-corrected H α , and CO (2–1) emission at 70–160 pc resolution in the first four PHANGS–JWST targets. We report correlation strengths, intensity ratios, and power-law fits relating emission in JWST’s F770W, F1000W, F1130W, and F2100W bands to CO and H α . At these scales, CO and H α each correlate strongly with mid-IR emission, and these correlations are each stronger than the one relating CO to H α emission. This reflects that mid-IR emission simultaneously acts as a dust column density tracer, leading to a good match with the molecular-gas-tracing CO, and as a heating tracer, leading to a good match with the H α . By combining mid-IR, CO, and H α at scales where the overall correlation between cold gas and star formation begins to break down, we are able to separate these two effects. We model the mid-IR above I ν = 0.5 MJy sr −1 at F770W, a cut designed to select regions where the molecular gas dominates the interstellar medium (ISM) mass. This bright emission can be described to first order by a model that combines a CO-tracing component and an H α -tracing component. The best-fitting models imply that ∼50% of the mid-IR flux arises from molecular gas heated by the diffuse interstellar radiation field, with the remaining ∼50% associated with bright, dusty star-forming regions. We discuss differences between the F770W, F1000W, and F1130W bands and the continuum-dominated F2100W band and suggest next steps for using the mid-IR as an ISM tracer. 

We measure empirical relationships between the local star formation rate (SFR) and properties of the star-forming molecular gas on 1.5 kpc scales across 80 nearby galaxies. These relationships, commonly referred to as “star formation laws,” aim at predicting the local SFR surface density from various combinations of molecular gas surface density, galactic orbital time, molecular cloud free fall time, and the interstellar medium dynamical equilibrium pressure. Leveraging a multiwavelength database built for the Physics at High Angular Resolution in Nearby Galaxies (PHANGS) survey, we measure these quantities consistently across all galaxies and quantify systematic uncertainties stemming from choices of SFR calibrations and the CO-to-H₂conversion factors. The star formation laws we examine show 0.3–0.4 dex of intrinsic scatter, among which the molecular Kennicutt–Schmidt relation shows a ∼10% larger scatter than the other three. The slope of this relation rangesβ≈ 0.9–1.2, implying that the molecular gas depletion time remains roughly constant across the environments probed in our sample. The other relations have shallower slopes (β≈ 0.6–1.0), suggesting that the star formation efficiency per orbital time, the star formation efficiency per free fall time, and the pressure-to-SFR surface density ratio (i.e., the feedback yield) vary systematically with local molecular gas and SFR surface densities. Last but not least, the shapes of the star formation laws depend sensitively on methodological choices. Different choices of SFR calibrations can introduce systematic uncertainties of at least 10%–15% in the star formation law slopes and 0.15–0.25 dex in their normalization, while the CO-to-H₂conversion factors can additionally produce uncertainties of 20%–25% for the slope and 0.10–0.20 dex for the normalization.

 

Abstract We measure the low- J CO line ratios R 21 ≡ CO (2–1)/CO (1–0), R 32 ≡ CO (3–2)/CO (2–1), and R 31 ≡CO (3–2)/CO (1–0) using whole-disk CO maps of nearby galaxies. We draw CO (2–1) from PHANGS-ALMA, HERACLES, and follow-up IRAM surveys; CO (1–0) from COMING and the Nobeyama CO Atlas of Nearby Spiral Galaxies; and CO (3–2) from the James Clerk Maxwell Telescope Nearby Galaxy Legacy Survey and Atacama Pathfinder Experiment Large APEX Sub-Millimetre Array mapping. All together, this yields 76, 47, and 29 maps of R 21 , R 32 , and R 31 at 20″ ∼ 1.3 kpc resolution, covering 43, 34, and 20 galaxies. Disk galaxies with high stellar mass, log ( M ⋆ / M ⊙ ) = 10.25 – 11 , and star formation rate (SFR) = 1–5 M ⊙ yr −1 , dominate the sample. We find galaxy-integrated mean values and a 16%–84% range of R 21 = 0.65 (0.50–0.83), R 32 = 0.50 (0.23–0.59), and R 31 = 0.31 (0.20–0.42). We identify weak trends relating galaxy-integrated line ratios to properties expected to correlate with excitation, including SFR/ M ⋆ and SFR/ L CO . Within galaxies, we measure central enhancements with respect to the galaxy-averaged value of ∼ 0.18 − 0.14 + 0.09 dex for R 21 , 0.27 − 0.15 + 0.13 dex for R 31 , and 0.08 − 0.09 + 0.11 dex for R 32 . All three line ratios anticorrelate with galactocentric radius and positively correlate with the local SFR surface density and specific SFR, and we provide approximate fits to these relations. The observed ratios can be reasonably reproduced by models with low temperature, moderate opacity, and moderate densities, in good agreement with expectations for the cold interstellar medium. Because the line ratios are expected to anticorrelate with the CO (1–0)-to-H 2 conversion factor, α CO 1 − 0 , these results have general implications for the interpretation of CO emission from galaxies. 

We combine JWST observations with Atacama Large Millimeter/submillimeter Array CO and Very Large Telescope MUSE Hαdata to examine off-spiral arm star formation in the face-on, grand-design spiral galaxy NGC 628. We focus on the northern spiral arm, around a galactocentric radius of 3–4 kpc, and study two spurs. These form an interesting contrast, as one is CO-rich and one CO-poor, and they have a maximum azimuthal offset in MIRI 21μm and MUSE Hαof around 40° (CO-rich) and 55° (CO-poor) from the spiral arm. The star formation rate is higher in the regions of the spurs near spiral arms, but the star formation efficiency appears relatively constant. Given the spiral pattern speed and rotation curve of this galaxy and assuming material exiting the arms undergoes purely circular motion, these offsets would be reached in 100–150 Myr, significantly longer than the 21μm and Hαstar formation timescales (both < 10 Myr). The invariance of the star formation efficiency in the spurs versus the spiral arms indicates massive star formation is not only triggered in spiral arms, and cannot simply occur in the arms and then drift away from the wave pattern. These early JWST results show that in situ star formation likely occurs in the spurs, and that the observed young stars are not simply the “leftovers” of stellar birth in the spiral arms. The excellent physical resolution and sensitivity that JWST can attain in nearby galaxies will well resolve individual star-forming regions and help us to better understand the earliest phases of star formation.

 

We present a high-resolution view of bubbles within the Phantom Galaxy (NGC 628), a nearby (∼10 Mpc), star-forming (∼2M_⊙yr⁻¹), face-on (i∼ 9°) grand-design spiral galaxy. With new data obtained as part of the Physics at High Angular resolution in Nearby GalaxieS (PHANGS)-JWST treasury program, we perform a detailed case study of two regions of interest, one of which contains the largest and most prominent bubble in the galaxy (the Phantom Void, over 1 kpc in diameter), and the other being a smaller region that may be the precursor to such a large bubble (the Precursor Phantom Void). When comparing to matched-resolution Hαobservations from the Hubble Space Telescope, we see that the ionized gas is brightest in the shells of both bubbles, and is coincident with the youngest (∼1 Myr) and most massive (∼10⁵M_⊙) stellar associations. We also find an older generation (∼20 Myr) of stellar associations is present within the bubble of the Phantom Void. From our kinematic analysis of the HI, H₂(CO), and Hiigas across the Phantom Void, we infer a high expansion speed of around 15 to 50 km s⁻¹. The large size and high expansion speed of the Phantom Void suggest that the driving mechanism is sustained stellar feedback due to multiple mechanisms, where early feedback first cleared a bubble (as we observe now in the Precursor Phantom Void), and since then supernovae have been exploding within the cavity and have accelerated the shell. Finally, comparison to simulations shows a striking resemblance to our JWST observations, and suggests that such large-scale, stellar-feedback-driven bubbles should be common within other galaxies.